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Abstract 

The Configuration Interaction (CI) method is applied to the calculation of the structures of a 
number of positron binding systems, including e + Be, e + Mg, e + Ca and e + Sr. These calculations 
were carried out in orbital spaces containing about 200 electron and 200 positron orbitals up to 
£ = 12. Despite the very large dimensions, the binding energy and annihilation rate converge 
slowly with £, and the final values do contain an appreciable correction obtained by extrapolating 
the calculation to the £ — > oo limit. The binding energies were 0.00317 hartree for e + Be, 0.0170 
hartree for e + Mg, 0.0189 hartree for e + Ca, and 0.0131 hartree for e + Sr. 

PACS numbers: 36.10.-k, 36.10.Dr, 31.25.Eb, 34.85. +x 



*Electronic address: mbromley@physics.sdsu.edu 
^Electronic address: jxml07@rsphysse.anu.edu.ai 



1 



I. INTRODUCTION 



The ability of positrons to bind to a number of atoms is now well established 1 , 2 , |3j], and 
all of the group II elements of the periodic table are expected to bind a positron l|, \$ . There 
have been two sets of calculations that are consistent, in that they tend to predict the same 
binding energy and annihilation rate. The first set of calculations were those undertaken 
on e + Be and e + Mg jl, H, Q with the fixed core stochastic variational method (FCSVM) 
PI El HI- Some time later, con figur ation interaction (CI) calculations were undertaken on 
e + Be, e + Mg, e + Ca and e + Sr (UtHl- The calculations for e + Be and e + Mg agreed to within 
the respective computational uncertainties, which were roughly about 5-10% for the binding 
energy. 

One feature common to all the CI calculations is the slow convergence of the binding 
energy and the annihilation rate. The attractive electron-positron interaction leads to the 
formation of a Ps cluster (Le. something akin to a positronium atom) in the outer valence 
region of the atom 0, 0, 0, The accurate representation of a Ps cluster using only 



single particle orbitals centered on the nucleus requires the inclusion of orbitals with much 
higher angular momenta than a roughly equivalent electron-only calculation 11, k| 3 3- 



For example, the largest CI calculations on the group II positronic atoms and PsH have 
typically have involved single particles bases with 8 radial function per angular momenta, £, 
and inclusion of angular momenta up to L max = 10 0,0, Even with such large orbital 



basis sets, between 5-60% of the binding energy and some 30-80% of the annihilation rate 
were obtained by extrapolating from L m „ = 10 to the L max = oo limit. 

Since our initial CI calculations [3, [lOj, [l6j , advances in computer hardware mean larger 
dimension CI calculations are possible. In addition, program improvements have removed 
the chief memory bottleneck that previously constrained the size of the calculation. As a 
result, it is now appropriate to revisit the group II atoms to obtain improved estimates of 
their positron binding energies and other expectation values. The new calculations that we 
have performed have orbital spaces more than twice as large as those reported previously. 
The estimated CI binding energies for all systems have increased, and furthermore the 
uncertainties resulting from the partial wave extrapolation have decreased. 

II. CALCULATION METHOD 

The CI method as applied to atomic systems with two valence electrons and a positron has 
been discussed previously jH^J, and only a brief description is given here. All calculations 
were done in the fixed core approximation. The effective Hamiltonian for the system with 
N e = 2 valence electrons and a positron was 



H = -\ v o - E \^ - ^ r °) + ^( r °) 

i=l 

+ V^dirW + KxcN + V p i{ri)) - V — 
i=i i=i 
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The index denotes the positron, while i and j denote the electrons. The direct potential 
(Vdir) represents the interaction with the electronic core, which was derived from a Hartree- 
Fock (HF) wave function of the neutral atom ground state. The exchange potential (V exc ) 
between the valence electrons and the HF core was computed without approximation. 

The one-body and two-body polarization potentials (V p i and V v <i) are semi-empirical with 
the short-range cut-off parameters derived by fitting to the spectra of their singly ionized 
ions. All details of the core-polarization potentials including the polarizabilities, a^, are 
given in Note that the functional form of the polarization potential, V p x, was set to 

be the same for the electrons and the positron. 

The positronic atom wave function is a linear combination of states created by multiplying 
atomic states to single particle positron states with the usual Clebsch-Gordan coupling 
coefficients ; 

\*;LS) = ^ Cij - (LiM^mjlLM^iSiMs^HjlSMs) 

x $i(Atom; LiSityfo) . (2) 

In this expression $>i(Atom; LiSi) is an antisymmetric atomic wave function with good L 
and S quantum numbers. The function <pj (r ) is a single positron orbital. The single particle 
orbitals are written as a product of a radial function and a spherical harmonic: 

0(r) = P(r)Y lm (r) . (3) 

As the calculations were conducted in a fixed core model we used HF calculations of the 
neutral atom ground states to construct the core orbitals. These HF orbitals were computed 
with a program that can represent the radial wave functions as a linear combination of Slater 
Type Orbitals (STO) 0. 

A linear combination of STOs and Laguerre Types Orbitals (LTOs) was used to describe 
the radial dependence of electrons occupying orbitals with the same angular momentum 
as those in the ground state. Orbitals that did not have any core orbitals with the same 
angular momentum were represented by a LTO set with a common exponential parameter. 
The STOs give a good representation of the wave function in the interior region while the 
LTOs largely span the valence region. The LTO basis 0, [Hj has the property that the 
basis can be expanded toward completeness without introducing any linear independence 
problems. 

The CI basis included all the possible L = configurations that could be formed by 
letting the two electrons and positron populate the single particle orbitals subject to two 
selection rules, 

max(4,4,£ 2 ) < L max , (4) 
mm(£ 1} £ 2 ) < L int . (5) 

In these rules £q is the positron orbital angular momentum, while l\ and ii are the angular 
momenta of the electrons. A large value of L max is necessary as the attractive electron- 
positron interaction causes a pileup of electron density in the vicinity of the positron. The 
Li n t parameter was used to eliminate configurations involving the simultaneous excitation 
of both electrons into high i states. Calculations on PsH and e + Be had shown that the 
choice of L; nt = 3 could reduce the dimension of the CI basis by a factor of 2 while having 
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an effect of about 1% upon the binding energy and annihilation rate |9(. The present set of 
calculations were all performed with L- in t = 4. 

Various expectation values were computed to provide information about the structure of 
these systems. The mean distance of the electron and positron from the nucleus are denoted 
by (r e ) and (r p ). The 27 annihilation rate for annihilation with the core and valence electrons 
was computed with the usual expressions [HI, 0, The 27 rate for the core (T c ) and 
valence (T,,) electrons are tabulated separately. 



A. Extrapolation issues 



The feature that differentiates mixed electron-positron CI calculations from purely elec- 
tron CI calculations is the slow convergence of the calculation with respect to L max , the 
maximum i of any electron or positron orbital included in the CI basis. Typically, a cal- 
culation is made to L max « 10 (or greater), with various extrapolation techniques used to 
estimate the L max — > 00 correction. For any expectation value one can write formally 

(X) Lmax = AxL > ( 6 ) 

L=0 

where AX L is the increment to the observable that occurs when the maximum orbital 
angular momentum is increased from L— 1 to L, e.g. 

AX L = (X) L - {X) L - 1 . (7) 

Hence, one can write formally 

00 

(x) = (x) Lmax + AxL ■ ( 8 ) 

L= L ma x4"l 

However, it is quite easy to make substantial errors in estimating the L max — > 00 correction 
There have been a number of investigations of the convergence of CI expansions 

for electronic and mixed electron-positron systems Q, H3, E3, HE ISS E3, ISE H, H3, E| . The 

reliability of the different methods to estimate the L max — > oo correction for the energy and 
annihilation rate has been assessed in detail elsewhere 2^1 . In this work, only the briefest 
description of the recommended methods are described. 



The recent computational investigations of helium [3JJ] and some positron-atom systems 
22|| suggest that usage of an inverse power series of the generic type 



(-^max 5) (-^max 



+ ~r DX ^ + , +--- > 0) 



(■^max 

is the best way to determine the L max — > 00 correction for the energy, E and the 27 annihi- 
lation rate. A three term series with n = 4 is used for the energy. One needs four successive 
values of E L to determine the coefficients Be, Ce and De- Once the coefficients have been 
fixed, the inverse power series is summed to J max = 100, after which the approximate result 

00 1 1 

L + (p-lX^ax + l)^ 1 ' (10) 
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is used |23j ]. 

The correction to T follows the same general procedure as the energy, but with two 
differences. The power in eq. © is set to n = 2 and only 2-terms are retained in the series 
(requiring three successive values of T L ). 

The usage of the inverse power series is the preferred approach when the asymptotic form 
for AX L has been established by perturbation theory. For other operators it is best to to 
use a single-term inverse power series with an indeterminate power, e.g 

AX L = - A - . (11) 
The factors A and p can be determined from the three largest calculations using 

'-(^)/ ta fe3)' 

and 

A = AX imax (L max + i) p . (13) 

Once p and A are determined, the L max — > oo correction can be included using the same 
procedure as adopted for the multi-term fits to the energy and annihilation. This method is 
used in determination of the L max — > oo estimates of (r e ), (r p ) and T c . However, the value 
of p is computed for all operators since it is useful to know whether pe and pr v are close to 
the expected values of 4 and 2 respectively. While the subdivision of the annihilation rate 
into core and valence components is convenient for physical interpretation, it was also done 
on mathematical grounds. The calculation of T c does not explicitly include correlations 
between the core electrons and the positron, and so the AI^ increments converge faster 
than the AI^ increments (i.e. pr c > Pr v )- 



III. CALCULATION RESULTS 

A. Improved FCSVM data for e+Be and e + Mg 

The FCSVM jj], |(| has also been applied to determine the structures of e + Be and e + Mg 
0,0. The FCSVM expands the wave function as a linear combination of explicitly correlated 
gaussians (ECGs), with the core orbitals taken from a HF calculation. One- and two-body 
polarization potentials are included while orthogonality of the active electrons with the core 
is enforced by the use of an orthogonalizing pseudo-potential 0, H, H3 • The FCSVM model 
hamiltonians are very similar to those used in the CI calculations. But there are some small 
differences in detail that lead to the FCSVM hamiltonian giving slightly different energies. 

The best previous FCSVM wave function for e + Be 0| gave a binding energy, 0.03147 
hartree, and annihilation rate 0.420 x 10 9 sec -1 , that were close to convergence. Some exten- 
sive re-optimizations seeking to improve the quality of the wave function in the asymptotic 
region yielded only minor changes (of the order of 1%) in the ground state properties [33]. 
Nevertheless, the latest energies and expectation values for the e + Be ground state are tab- 
ulated in Tables and HIl These values should be converged to better than 1% with respect 
to further enlargement and optimization of the ECG basis. 

The more complex core for Mg does slow the convergence of the energy and other prop- 
erties of e + Mg considerably [l|. The best energy previously reported for this system was 
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0.016096 hartree |33j|. The current best wave function, which is constructed from a linear 
combination of 1200 ECGs gives a binding energy of 0.016930 hartree and a valence annihila- 
tion rate of 1.0137 x 10 9 sec^ 1 . Other expectation values are listed in Table UJ Examination 
of the convergence pattern during the series of basis set enlargements and optimizations 
suggests that the binding energy and annihilation rate are converged to between 2% and 
5%. 

The FCSVM binding energies do have a weak dependence on one parameter in the cal- 
culation since the orthogonalizing pseudo-potential is actually a penalty function, viz 

\P= MhKM' ( 14 ) 

i £ core 

that was added to the hamiltonian. Choosing A to be large and positive means the energy 
minimization automatically acts to construct a wave function which has very small overlap 
with the core jy, |32j, |34| . The FCSVM properties reported in Tables HI and ITT1 were computed 
with A = 10 5 hartree. The core overlap (i.e. the expectation value of P) was 1.86 x 10~ n 
for e+Be and 1.61 x 10~ 10 for e+Mg. 

B. CI results for group II atoms 

Table H] contains the results of the current series of calculations on the four positronic 
atoms. The size of the calculations for the four atoms were almost the same. The electron- 
electron angular momentum selector was set to L int = 4. For i > 3 at least 15 LTOs 
were included in the radial basis sets for the electron and positron orbitals. For £ < 2 the 
dimension of the orbital basis sets were slightly larger than 15 and the basis sets for electrons 
occupying orbitals with the same angular momentum as those in the core were typically a 
mix of STOs (to describe the electron close to nucleus) and LTOs. The calculations used 
basis sets with L max = 9, 10, 11 and 12. The calculations with L max < 12 had configuration 
spaces which were subsets of the L max = 12 and this expedited the computations since one 
list of radial matrix elements was initially generated for the L max = 12 basis and then reused 
for the smaller basis sets. 

The secular equations that arose typically had dimensions of about 500,000 and the 
diagonalizations were performed with the Davidson algorithm using a modified version of 
the program of Stathopolous and Froese-Fischer j3|| . Convergence was not very quick and 
about 16000 iterations were needed to achieve convergence in some cases. It was possible 
to speed up the diagonalization for L max < 12. An edited eigenvector from the L max = 12 
calculation was used as the initial eigenvector estimate, and this often reduced the number 
of iterations required by 50%. 

1. Results for e + Be 

The lowest energy dissociation channel is the e + + Be channel, which has an energy of 
— 1.01181167 hartree with respect to the doubly ionized Be 2+ core. The agreement of the 
extrapolated CI binding energy of e — 0.003169 hartree with the FCSVM binding energy 
of e — 0.003180 is better than 1%. A similar level of agreement exists for the (r e ) and (r p ) 
expectation values. 
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The only expectation value for which 1% level of agreement does not occur is the anni- 
hilation rate and here the extrapolated CI value of 0.4110 x 10 9 sec -1 is only about 3.5% 
smaller than the FCSVM of 0.4267 x 10 9 sec -1 . However, it is known that the convergence of 
the annihilation rate with respect to an increasing number of radial basis functions is slower 
than the convergence of the energy [i^, HDJ • This means that a CI type calculation has an 
inherent tendency to underestimate the annihilation rate. For example, a CI calculation on 
PsH of similar size to the present e + Be calculation underestimated the annihilation rate by 
6% 2^1 . That the exponent of the polar law decay, pr v = 2.10, is larger than the expected 



asymptotic value of p = 2.0 is consistent with this idea. A better estimate of the annihilation 
rate can be obtained by simply forcing Cr to be zero in eq. and thus using Ar 12 to fit 
-Br- When this is done done the annihilation rate increases to 0.4178 x 10 9 sec -1 . 



2. Results for e + Mg 

The results of the calculations with e + Mg are listed in Table |U The lowest energy 
dissociation channel is to e + + Mg, which has an energy of —0.83291427 hartree with respect 
to the doubly ionized Mg 2+ core. 

The CI calculations, reported in Table 0for L max = 9, 10, 11 and 12 are largely consistent 
with the FCSVM calculations. The largest explicit CI calculation gives a binding energy 
of 0.015658 hartree. Extrapolation to the L max — > oo limit adds about 10% to the binding 
energy, and the final estimate was 0.017040 hartree. Despite the better than 1% agreement 
between the CI and FCSVM calculations, a further binding energy increase of about 1-2% 
would be conceivable if both calculations were taken to the variational limit. 

The slow convergence of Y v with L max is evident from Table U and the extrapolation 
correction contributes about 36% to the overall annihilation rate. The present L max — > oo 
estimate can be expected to be too small by 5-10%. 

All the other expectation values listed in Table lie with 1-2% of those of the FCSVM 
expectation values. As a general rule, inclusion of the L max — > oo corrections generally 
improves the agreement between the CI and FCSVM calculations. 



3. Results for e + Ca 

The results of the calculations with e + Ca are listed in Table HI Since neutral calcium 
has an ionization potential smaller than the energy of Ps ground state (the present model 
potential and electron orbital basis gives -0.43628656 hartree for the Ca + energy and - 
0.65966723 hartree for the neutral Ca energy), its lowest energy dissociation channel is the 
Ps + Ca + channel. The present model potential gives this channel an energy of —0.68628656 
hartree. 

The energies listed in Table H] indicate that e + Ca is the positronic atom with the largest 
known binding energy, namely e = 0.018929 hartree. The L max — > oo correction contributes 
20 % of the binding energy. The partial wave series is more slowly convergent for e + Ca than 
for e + Mg (i.e. pe is smaller, and the coefficients Ce and D E in eq. (jTUj) are larger). This is 
expected since calcium has a smaller ionization potential, and so the electrons are located a 
greater distance away from the nucleus. This makes it easier for the positron to attract the 
electrons, and the stronger pileup of electron density around the positron further from the 
nucleus requires a longer partial wave expansion to represent correctly. 
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The slower convergence of the wave function with L max makes an even larger impact on 
the annihilation rate. Some 41% of the annihilation rate of T v = 1.478 x 10 9 sec -1 comes 
from the L max — > oo correction. As mentioned earlier for e + Mg, it is likely that this value is 
slightly smaller than the true annihilation rate. 

The extrapolation corrections for (r p ) and r c listed in Table are unreliable. The e + Ca 
system, at large distances consists of Ca + + Ps. In other calculations of positron binding 
systems it has been noticed that systems that decay asymptotically into Ps + X do not have 
an (r p ) that changes monotonically with L max |9j, [16J . Initially, the positron becomes more 
tightly bound to the system as L max increases, resulting in a decrease in (r p ). However, (r p ) 
tends to increase at the largest values of L max . The net result of all this is that A(r p ) L (and 
by implication Ar^) approach their asymptotic forms very slowly. The best policy is to 
simply not to give any credence to the extrapolation corrections for either of these operators 
for e + Ca (and e + Sr). The small value of p for A(r e ) L suggests that the reliability of the 
L mm — > oo correction may be degraded for this expectation value as well. 

4- Results for e + Sr 

The results of the calculations for e + Sr are listed in Table|lJ Since neutral strontium has an 
ionization potential smaller than the energy of Ps ground state (the present model potential 
and electron orbital basis gives -0.40535001 hartree for the Sr + energy and -0.61299101 
hartree for the neutral Sr energy), its lowest energy dissociation channel is the Ps + Sr + 
channel, which has an energy of -0.65535001 hartree. The small ionization potential of 
0.20764100 hartree means that the structure of the e + Sr ground state will be dominated by 
a \I/(Sr + )\l/(Ps) type configuration [l[. This leads to slower convergence of the ground state 
with L max which is evident from Table |U 

As expected, the binding energy of e + Sr is smaller than that of e + Ca. Previous inves- 
tigations have indicated that positron binding energies should be largest for atoms with 
ionization potentials closest to 0.250 hartree (the Ps binding energy) 4, 36]. There is obvi- 
ously some uncertainty in the precise determination of the binding energy due to fact that 
I max — > oo correction constitutes some 37% of the binding energy of 0.013102 hartree. The 
net effect of errors due to the extrapolation correction are not expected to be excessive. Ap- 
plying eq. © with only the first two-terms retained (i.e. De = 0) results in a final energy 
0.012764 hartree, which is 3% smaller than the value of 0.013102 Hartree. The present e + Sr 
binding energy is some 30% larger than the energy of the previous CI calculation listed in 
Table 

The final estimate of the valence annihilation rate was 1.553 x 10 9 sec -1 and some 43% 
of the annihilation rate comes from the L max — > oo correction. This value of T v could easily 
be 10% smaller than the true annihilation rate. The explicitly calculated expectation values 
for (r e ), (r p ) and T c at L max = 12 should be preferred since the L max — > oo corrections in 
these cases are likely to be unreliable. 

C. 3-body clustering 

While the truncation of the basis to L int = 4 has little effect on the e + Be system, its 
effect is larger for the e + Sr system. The more loosely bound alkaline-earth atoms have their 
electrons localized further away from the nucleus, and this makes it easier for the positron 
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to form something like a Ps~ cluster |4j, |37( • When this occurs, correlations of the positron 
with both electrons increase in strength, and the inclusion of configurations with L mt > 4 
becomes more important. 

The relative size of of these neglected L int > 4 configurations can be estimated using 
techniques similar to those adopted for the L max — > oo corrections. Calculations for a 
succession of L int values were performed in earlier works 0, Hoi ]. The assumption is made 
that the binding energy and annihilation rate increments scale as A/(L int + |) 4 (note, the 
power of 4 for the annihilation is used since L- mt only has a directly effect on electron-electron 
correlations). The difference between an L int = 2 and L int = 3 is used to estimate A and 
then eq.( ITU]) determines the L int — > oo correction (in the case of e + Be calculations up to 
L int = 10 exist @). 

Table |H] contains a summary of the final binding energies obtained from the present CI 
calculations, and earlier binding energies obtained alternate methods. As part of this table, 
energies with an additional L int — > oo correction are also given. The size of the correction 
ranges from 1.8 x 10~ 5 hartree for e + Be to 21.9 x 10~ 5 hartree for e + Sr. Even though these 
estimations of the correction are not rigorous, they indicate that the underestimation in the 
binding energy resulting from a truncation of the configuration space to L int < 4 is most 
likely to be 2% or smaller. 

A similar analysis could be done for the annihilation rate but previous results indicate 
that r„ is less sensitive than e to an increase in Lj nt 0,^1 . The net increases in Y v for e + Be, 
e+Mg, e+Ca and e+Sr were 0.0011 x 10 9 sec" 1 , 0.0030 x 10 9 sec" 1 , 0.0039 x 10 9 sec" 1 and 
0.0039 x 10 9 sec -1 respectively. All of these extra contributions to Y v correspond to changes 
of less than 0.5%. 



IV. SUMMARY AND CONCLUSIONS 

The summary of binding energies, produced by the current methods and other completely 
different approaches presented in Table |H] shows that the only methods that consistently 
agree with each other are the CI and FCSVM calculations. Both these methods are varia- 
tional in nature, both use realistic model potentials designed on very similar lines, and both 
have shown a tendency for the bin ding energies to slowly creep upwards as the calculation 
size is increased (refer to refs |l|, |y, [lOj for examples of earlier and slightly smaller binding 
energies). The PO and MBPT approaches do not give reliable binding energies. 

The diffusion Monte Carlo method |38| gives an e + Mg binding energy of 0.0168±0.0014 
hartree which is very close to the present energy. This calculation was fully ab-initio and 
did not use the fixed core approximation. However, application of the same diffusion Monte 
Carlo method to e + Be gave a binding energy which is only half the size of the present value. 

The present binding energies are all larger than those given previously 0, due to 
the usage of a radial basis which was almost twice the size of earlier calculations. In two 
cases, + Ca and e + Sr the increase in binding energy exceeds 10%. The binding energies for 
e + Be and e + Mg are in agreement with those of FCSVM calculations to within their mutual 
uncertainties. Further enlargement of the basis could lead to the positron binding energies 
for Mg, Ca and Sr increasing by a few percent. 

Estimates of the annihilation rate have also been extracted from the CI wave functions. 
The present annihilation rates are certainly underestimates of the true annihilation rate. 
The annihilation rate converges very slowly with respect to the radial basis and similar 
sized calculations on PsH suggest that the present annihilation rates could easily be too 
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small by at least 5% HJ E2 . 

The speed at which the partial wave expansion converges with respect to L max is seen 
to decrease as the ionization energy of the parent atom decreases jl|, • In addition, the 
importance of 3-body clustering (i.e. convergence with respect to L int ) was seen to increase 
as the ionization energy of the parent atom decreased 0. 

The main factor limiting the size of the calculations now is the time taken to perform 
the diagonalizations. Although, the calculations were performed on a Linux/Myrinet-based 
cluster, the sheer number of iterations, (16000 in the worst case), used by the Davidson 
method, meant that it could take 30 days to perform a diagonalization using 24 CPUs. 
However, the main reason for a dop ting the Davidson method was the availability of a pro- 
gram that was easy to modify [35| . Usage of the more general Lanczos method |43[ might 
lead to a quicker diagonalization and thus permit even larger calculations. 
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TABLE I: Results of CI calculations for positronic alkaline-earth atoms for a given -L max . The 
E column gives the three-body energy with respect to the doubly ionized frozen core and e is the 
binding energy with respect to the lowest energy dissociation channel. The T v and r c columns give 
the valence and core annihilation rate (in 10 9 sec -1 ) . The results in the row 10* are taken from 
earlier CI calculations for these systems with L max = 10 (and L- in t = 3). The results in the 

row oo use the methods described in the body of the text to evaluate the L max — > oo correction. 
The exponent p characterizes the rate of decay of the expectation value increments evaluated at 
£max = 12 using eq. (JT2J). 



T 

-^max 




£ 


(re) 


( r p) 


T"! 

r c 










e + Be 








10* 

9 


(-1.0143769) 
-l. 01435756 


(0.002533) 
0.00254589 


(2.639) 
2.6388477 


(10.746) 
10.874256 


(0.001962) 
0.00193993 


(0.2411) 
0.24026720 


10 


-I f~\ t A A o -I r^i 

-1.01448318 


0.00267151 


2.6418168 


10.699433 


0.00198405 


0.25651443 


11 


-1.01457837 


0.00276670 


2.6441227 


10.574208 


0.00201619 


0.27004634 


12 


-1.01465138 


0.00283971 


2.6459282 


10.482126 


0.00204005 


0.28140404 


P 


3.1806 


3.1806 


2.9339 


3.6871 


3.5764 


2.1006 


oo 


-l. 0149809 


0.0031692 


2.65673 


10.09755 


c\ nnoi A A 

0.002144 


0.410976 


FCSVM 


-1.0151335 


0.003180 


2.654 


10.048 


0.00221 


0.4267 


10* 
9 


(-0.8473592) 

/~\ o A r—? A -i A r\ A 

-0. 84741494 


(0.0145092) 
0. 01450067 


e + Mg 
(3.382) 
3.3831320 


(7.101) 

1 1 /Tin 

7.116532 


(0.010845) 
0. 01079647 


(0.5429) 

r\ r < AAAA1 A 

0. 54089010 


10 


-0.84790548 


0.01499121 


3.3936654 


7.071950 


f\ f\~\ f\Ci Af\ A A 

0.01084944 


0.57692369 


11 


-0.84828090 


0.01536663 


3.4022694 


7.040929 


0.01087921 


0.60775407 


12 


-0.84857204 


0.01565777 


3.4093312 


7.018703 


0.01089568 


0.63435278 


P 


3.0496 


3.0496 


, -j r* r\r\ 

2.3690 


3.9985 


7.0953 


1.7706 


oo 


-0.8499543 


0.0170400 


3.47039 


6.93657 


0.010922 


0.990069 


FCSVM 


-0.849002 


0.016930 


O A A T 

3.447 


6.923 


0.0112 


1.0137 


10* 
9 


(-0.6986443) 
-0.69855551 


(0.0123578) 
0.01226895 


e + Ca 
(4.456) 
4.4602428 


(6.848) 
6.863740 


(0.01355) 
0.01343426 


(0.7335) 
0.72709017 


10 


-0.69975764 


0.01347109 


4.4873848 


6.872414 


0.01323075 


0.78001274 


11 


-0.70069553 


0.01440898 


4.5110869 


6.885039 


0.01304512 


0.82640757 


12 


-0.70143637 


0.01514981 


4.5315631 


6.898804 


0.01288316 


0.86733542 


P 


2.8286 


2.8286 


1.7546 


-1.0371 


1.6361 


1.5037 


oo 


-0.7052160 


0.0189295 


4.86076 




0.009780 


1.478148 


10* 
9 


(-0.6602186) 
-0.65997599 


(0.0048689) 
0.00462598 


e + Sr 
(4.850) 
4.8638673 


(7.056) 
7.100141 


(0.01487) 
0.01464684 


(0.7488) 
0.73239378 


10 


-0.66146709 


0.00611708 


4.8979559 


7.123685 


0.01432317 


0.78845209 


11 


-0.66263875 


0.00728874 


4.9283728 


7.150071 


0.01403253 


0.83790890 


12 


-0.66357065 


0.00822064 


4.9552753 


7.176708 


0.01377785 


0.88177286 


P 


2.7459 


2.7459 


1.4725 


-0.1134 


1.5844 


1.4393 


oo 


-0.6684520 


0.0131020 


5.65380 
12 




0.008456 


1.552589 



TABLE II: Binding energies (in hartree) of positronic beryllium, magnesium, calcium and stron- 
tium. Only the latest calculations of a given type by a particular group are listed in this table. 



Calculation 


e+Be 


+eMg 


e+Ca 


CI (L max = 12) 


0.002840 


0.015658 


0.015150 


CI (L max -> oo) 


0.003169 


0.017040 


0.018929 


CI (L int OO) 


0.003187 


0.017099 


0.019122 


Previous-CI a 


0.003083 


0.01615 


0.01650 


FCSVM 


0.003161 


0.016930 




DMC b 


0.0012(4) 


0.0168(14) 




SVM c 


0.001687 






PO d 




0.00055 




PO e 




0.00459 




MBPT / 




0.0362 





e + Sr 
0.008221 
0.013102 
0.013321 
0.01005 



"Previous CI (L max -> oo) [Ell 

h DMC, the statistical .uncertainty in the last digit(s) is given in the brackets [3 
c Fully ab-initio SVM [g 

^Polarized orbital calculation, jiipole only |39| 
e Polarized orbital calculation 40] 



■''Many Body perturbation theory [4l| 
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